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We discover that the multiple degrees of freedom associated with magnetic skyrmions: size, position, and he- 
licity, can all be used to control the Josephson effect and O-tt transitions occurring in superconductor/magnetic 
skyrmion/superconductor junctions. In the presence of two skyrmions, the Josephson effect depends strongly on 
their relative helicity and leads to the possibility of a helicity-transistor effect for the supercurrent where the crit¬ 
ical current is changed by several orders of magnitude simply by reversing the helicity of a magnetic skyrmion. 
Moreover, we demonstrate that the Fraunhofer pattern can show a local minimum at zero flux as a direct result 
of the skyrmion magnetic texture. These findings demonstrate the rich physics that emerges when combin¬ 
ing topological magnetic objects with superconductors and could lead to new perspectives in superconducting 
spintronics. 
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The interplay between superconductivity and ferromag¬ 
netism in hybrid structures has received much attention in 
recent yearsij^, due to its allure from a fundamental physics 
viewpoint and also because of improved and new functionality 
brought about by using superconductors in spintronics^. Due 
to the proximity effect, the Cooper pairs induced in the ferro- 
magnet acquire a finite center of mass momentum. There¬ 
fore, the pair amplitude oscillates in space which may re¬ 
sult in a sign change of the Josephson current in ferromag¬ 
netic Josephson junctions. This effect can be used to con¬ 
trol the quantum ground state of the system, altering from 
a state where the superconducting phase difference is 0 to a 
state where it is tt——. This O-tt transition was originally ob¬ 
served in Josephson junctions through weak ferromagnetsi^ 
Also, in the presence of inhomogeneity of the magnetic 
order, triplet pairing with spin aligned with the local ex¬ 
change field is generated in the ferromagnet due to spin 
flip scatterin g**^!** . Experiments have successfully demon¬ 
strated the presence of such spin-triplet pairing by observing 
a Josephson current through strong ferromagnetsi^— , which 
can be explained via the concepts of spin-mixing and spin- 
rotation taking place near the superconductor/ferromagnet 
interfacei^. Using equal spin triplet pairings, the possibil¬ 
ity arises to enhance existing effects or discover new ones in 
spintronicsii^^. Recently, it has been also proposed that in¬ 
homogeneous ferromagnet/superconductor junctions can cre¬ 
ate topological superconductivityi^TS^ 

Currently, much interest is garnered by magnetic skyrmions 
in chiral magnets^—. Such objects are characterized by 
a topologically protected spin configuration. Due to their 
peculiar magnetic structure, several intriguing phenomena 
have been discovered such as topological and skyrmion Hall 
effects22r22 and current-driven motion of skyrmion with ul- 
tralow current density^—. It has been shown that magnetic 
skyrmions can be also driven by a temperature gradient^^^. 
A thermal gradient is predicted to induce a skyrmion mo¬ 
tion towards the high temperature region accompanied by a 
skyrmion Hall effect^. Skyrmions are accompanied by a 
degree of freedom known as their helicity, which is deter¬ 
mined by their spin swirling direction. It has been experi¬ 


mentally demonstrated that the helicity of skyrmions can be 
changed both via a small external magnetic field^ and spin- 
orbit interactions^. This opens the exciting prospect that any 
physical quantity that responds to a change in the skyrmion 
helicity degree of freedom will be controllable via an external 
held. 

In this paper, we investigate how the degrees of freedom 
of magnetic skyrmions, such as helicity, can influence the 
supercurrent-response and quantum ground state of Joseph¬ 
son junctions including skyrmions. We find that the supercur¬ 
rent is strongly influenced by the (i) size, (ii) position, and (iii) 
helicity of magnetic skyrmions. We discover that the O-tt tran¬ 
sition can in fact be triggered by changing any of these three 
skyrmion properties, which in turn have been confirmed to be 
experimentally tunable via external magnetic fields^, spin- 
orbit interactions^, and electric currents^Ml. This offers a 
new and dynamical way of manipulating the quantum ground 
state of a superconducting system via magnetic skyrmions. 
We then show that the strong dependence on the helicity cre¬ 
ates a helicity-transistor effect for supercurrents, where the 
critical current is changed by several orders of magnitude 
upon reversing the helicity of a skyrmion. Moreover, we And 
that as a direct consequence of the skyrmion magnetization 
texture, the Fraunhofer pattern can display a local minimum at 
zero flux in contrast to conventional homogeneous magnetic 
Josephson junctions. In what follows, we will demonstrate 
these properties in Josephson junctions featuring both single 
and two skyrmions. 

We consider a 2D superconductor/ magnetic skyrmion / su¬ 
perconductor junction as shown in Fig. [T] This geometry is 
also expected to approximate well a planar junction geometry 
with two separated superconducting electrodes deposited on 
top of a thin magnetic film containing skyrmions. By assum¬ 
ing that the proximity effect is weak, we utilize the linearized 
Usadel equationji^dl 

DVV. - 2u;„/, - 2/ft • h = 0, (1) 

DVHt - 2w„ft - 2//,h = 0. (2) 

Here, D and a;„ are the diffusion constant in the magnet and 
Matsubara frequency, respectively, fs is the singlet anoma- 
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lous Green’s function while f( represents the triplet anoma¬ 
lous Green’s functions, h is the exchange field representing a 
magnetic structure with two skyrmions:— 


h = ^ 2 [2Re(u)x -f 2Im(u)y -f (1 - \u\'^)z], 

1 -f \u\ 

_ iX iX' 

x-Xc-i{y-yc) x - x'^ - i[y - y'^)' 


(3) 


Here, {xc,yc) and {x'^,y'^) determine the centers of the two 
skyrmions. A and A' are the characteristic sizes of the 
skyrmions. The signs of A and A' determine the helicities 
of the skyrmions. h is the magnitude of the exchange field. 
By setting A' = 0, the above exchange field represents a 
single skyrmion texture. We consider the magnetic region 
in —L/2 < x,y < Lj^. The interfaces are located at 
X = ±Lj2. 

The boundary condition at x = —L12 reads^ 

-lB^^ + Gsfs=Fs, + =0 (4) 

ox ox 

where fi (i = 1 , 2 ,3) is the components of it, and Gg and Fg 
are bulk Green’s functions in the superconductor given by 


Gg 


UJn 


A exp {—iipj2) 
+ A2 


(5) 


Here, 7 ^ describes the interface barrier strength, ^ is the su¬ 
perconducting coherence length, A is the gap function, and 
kp is the phase difference between the superconductors. The 
boudary condition at x = L /2 is given by changing the 
signs of the derivative and p in the above boundary condi¬ 
tion. The boundary condition at ?/ = ±L/2 reads = 0 
with a = s, 1, 2, 3. The Josephson current is calculated as 


elxF 

27rTc 


^ n>0 


(6) 


with the (transition) temperature T{Tc) and resistance of the 
magnet per length R. We define the total current as Ix = 
i Ixdy Sit X = — L/2. The critical current and that in¬ 

cluding the sign of the current are denoted by Ixc and I'xc’ 
respectively; Ixc = \I'xc\- B^low, we fix the parameters as 
7 b = 10 , T/Tc = 0.9 and /i/Aq = 1.5 where Ag denotes 
the gap energy at zero temperature. Calculation of the Joseph¬ 
son current requires a solution of the 2D Usadel equation. We 
have solved the Usadel equation numerically by using an iter¬ 
ative method. 

We begin by considering junctions with a single skyrmion 
(A' = 0) as shown in Fig. [Tla). In Fig. |2ja), we show the criti¬ 
cal current as a function of the length of the magnetic region L 
for several sizes of the skyrmion A. The skyrmion is assumed 
to be positioned in the center of the junction, Xc = yc = 0 . 
We find a O-tt transition as a function of L. It is also seen 
that the transition point can be controlled by altering the size 
of the skyrmion, A. In Fig. |2|b), we show the critical cur¬ 
rent as a function of A for several L and Xc = yc = 0. For 
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FIG. 1: (Color online) Schematic illustration of the superconduc¬ 
tor/magnet/superconductor junction with (a) a single skyrmion and 
(b) two skyrmions. This setup may also be viewed as a simpli¬ 
fied model for a lateral junction with superconducting electrodes de¬ 
posited on top of a skyrmion thin film. 



FIG. 2: (Color online) (a) The critical current as a function of the 
length of the magnetic region L for several sizes of the skyrmion 
A. (b) The critical current as a function of A for several L. We set 
Xc = yc = 0 and A' = 0. 
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FIG. 3: (Color online) The critical current including the sign of the 
current ^ function of the position of the skyrmion Xc and 

Uc for A/^ = 0.5, X' = 0, and L/^ = 4.2. This shows where the 
skyrmion should be located in the junction to induce a 0 —tt transition 
(the green circle). 



L/^ 



FIG. 4: (Color online) The critical current as a function of the mag¬ 
netic flux for several L. We set Xc = Vc = 0, A/^ = 0.5 and A' = 0. 


L/^ = 4, the O-TT transition occurs around A/^ = 0.42. These 
results indicate that the O-tt transition is tunable by chang¬ 
ing magnetic field and possibly also by applying an elec¬ 
tric field/gate voltage since this breaks inversion symmetry 
and hence can modify the Dzyaloshinskii-Moriya interaction 
which in turn changes the size of the skyrmion. The tunable 
size of skyrmions in helimagnetic alloys via spin-orbit cou¬ 
pling has been also experimentally verified in Ref4S, indicat¬ 
ing that the O-tt transition predicted here can be manipulated 
via changing the skyrmion size according this route. 

The next aspect we consider is how the skyrmion posi¬ 
tion influences the supercurrent response of the system. A 
unique feature of skyrmions is that the ultralow current den¬ 
sity (~ 10^ A/cm^) can induce their translational and/or ro¬ 
tational motions, which is typically 5 orders of magnitude 
smaller than the required density in conventional domain wall 
ferromagnets. This has been experimentally demonstrated in 
the helimagnet MnSi^ and FeGe^. Motivated by this, in Fig. 


FIG. 5: (Color online) The critical current as a function of the length 
of the magnetic region L for A/^ = 0.5 and either equal (A'/^ = 
0.5) or opposite (X '= —0.5) helicity of the two skyrmions. (a) 
Xc/^ = 0.5, xQ^ = —0.5 and yc = y'c = 0. (b) Xc/i = 1, 

x'a/l = -1 and i/c = ?/c = 0 . 


[3 we show the critical current including the sign of the cur¬ 
rent as a function of the position of the skyrmion Xc and yc 
for A/^ = 0.5 and L/^ = 4.2. It is found that a O-tt transi¬ 
tion occurs by changing the position of the skyrmion. Since 
the position of the skyrmion can be manipulated by current or 
temperature gradient, this offers a way to control the quantum 
ground state of the system. 

We also consider the critical current as a function of the 
magnetic flux threading through the magnetic region. We con- 
sidier uniform magnetic field along the z-axis and include the 
vector potential of the form A = B{—y, 0,0) by the substi¬ 
tution V —> V + i^Ain the equations. Figure |4] exhibits 
the critical current as a function of the magnetic flux $ for 
several L with $ = BL^ and $o = /i/2e. We find con¬ 
ventional Fraunhofer diffraction patterns for L/^ = 3.5 and 
5. However, for L/^ = 4.2 near the O-tt transition point as 
shown in Fig. 2(a), a minimum appears at <i> = 0 in the Fraun¬ 
hofer pattern. This can be understood as follows. Due to the 
skyrmion magnetization texture varying along the y-axis, the 
present junction may be regarded as a parallel circuit of 0 and 
TT junctions. In such a circuit of 0 and tt junctions, a local 
minimum at $ = 0 can appear due to the cancellation of the 
Josephson currents from the 0 and tt segments^. In this way, 
the Fraunhofer pattern in our setup can display a local mini¬ 
mum at zero flux near the O-tt transition points. 

Now, let us consider junctions with two skyrmions (see Fig. 
[nb)) and focus on the effect of the helicities. The presence of 
multiple skyrmions in the Josephson junction is particularly 
relevant in light of the experimental demonstration of multiple 
skyrmion configuration featuring skyrmions with both types 
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of helicities^2ii£. The helicity was shown to be reversible via 
an external field of order ~ 100 mT—. Figure|3shows the crit¬ 
ical current as a function of the length of the magnetic region 
L for A/^ = 0.5, considering both equal (A'/^ = 0.5) and op¬ 
posite (A'/f = —0.5) helicities of the two skyrmions. In Fig. 
|5](a), the positions of the skyrmions are set as Xc/^ = 0.5, 
x'^/^ = —0.5 and TJc = v'c = 0- Remarkably, we see that 
the O-TT transition point depends on the helicity, which means 
that the reversal of the helicity can induce a O-tt transition in 
itself. Moreover, this effect opens the possibility for a helicity- 
transistor effect for the supercurrent: close to a O-tt transition 
point for one helicity configuration, changing the relative he¬ 
licity in situ will result in an increase of the critical current 
density of several orders of magnitude as seen in Fig. E] In 
Fig. |5](b), the distance of the skyrmions is set to be longer 
as Xc/^ = 1, x'^/^ = — 1 and y^, = y'^ = 0. It is found 
that the O-tt transition points corresponding to the two helicity 
configurations become closer compared to Fig. |5](a). When 
two skyrmions are sufficiently separeted (e.g., for very large 
Xc), we can regard the two skyrmions as independent of each 
other and hence the effect of the helicity becomes negligible. 
Note that since the helicity changes under mirror operation 
with respect to xy plane, the results remain the same under 
the reversal of all the helicities of the skyrmions: the results 
depend on the relative sign of the helicities, and for a single 
skyrmion, the results do not depend on the helicity. 

For the chiral magnet MnSi, the material parameters are es¬ 
timated as h ^ 1 eV and A ~ 10nm22. The change of the 
exchange field h will shift the O-tt transiton point. Here, we 
have considered junctions with a single and two skyrmions. 
Skyrmions can also form a hexagonal lattice, and the appli¬ 


cation of our work to such a skyrmion configuration would 
be also informative. It could also be of interest to consider 
skyrmion tubes lined up along the junction direction, which 
should be experimentally feasible in layered thin-film struc¬ 
tures. Moreover, the unusually low threshold for current- 
induced skyrmion motion would be very interesting to investi¬ 
gate in the context of supercurrent-induced magnetization dy¬ 
namics and spin-transfer torques^. We leave these issues for 
future explorations. 

In summary, we have investigated the Josephson effect in 
superconductor/magnetic skyrmion/superconductor junction. 
It is found that the degrees of freedom associated with the 
skyrmions (size, position, and helicity), which recently have 
been demonstrated experimentally to be tunable via different 
routes, lead to a new dynamical way to control the O-tt tran¬ 
sition, offering the tantalizing prospect of a helicity-transistor 
for supercurrents. It is also shown that the Fraunhofer pattern 
can exhibit a local minimum at zero flux as a consequence of 
the skyrmion magnetization texture. 
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